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We develop an analytical formalism for studying optical analogues of spherically symmetric black- 
hole spacetimes. We demonstrate the exact similarity between the electromagnetic wave equations 
in an inhomogeneous medium in flat spacetime and those in a general-relativistic curved spacetime 
respectively. The permittivity and the permeability of an inhomogeneous optical medium act as 
metric components of the analogue spacetime. Manifest properties of black holes like curved light 
trajectories follow directly from our formalism. Other black-hole phenomena like quasinormal modes 
can also be studied within this framework. We extend the above formalism and obtain an approxi- 
mate analogue for the case of purely dielectric media. These optical black holes can be employed to 
investigate black-hole phenomena with table-top experiments. 

PACS numbers: 04.70.Bw,04.20.-q 
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I. INTRODUCTION: 

Advances in Transformation Optics and Analog Grav- 
ity have led to systems that closely resemble those en- 
gendered by the general theory of relativity such as the 
black holes. Models that have come under the category 
of Analog Gravity [lj make use of sound waves in flowing 
media Q like super fluids in Bose-Einstein Condensates 
[H [H to mimic certain aspects of a black hole . Such sys- 
tems developed until now are not exact analogues of any 
real spacetimes with corresponding metric tensors. On 
the other hand, transformation optics has utilized the 
formalism of general relativity Q to build gradient index 
static media in order to gain excellent control over the 
light paths in those media [6]. Although much of these 
efforts in transformation optics have been directed to- 
wards developing technological applications, the results 
can be used effectively to investigate and construct mod- 
els of optical analogues of gravity. In this context, we 
develop a general analytical formalism as the basis for 
optical black-hole analogues using static optical media. 
Our purpose is to generate in as much detail as possible, 
an optical analogue of the gravitational black hole and 
its manifestly characteristic phenomena, such as the exis- 
tence of a one-way membrane, namely the event horizon, 
the bending of light rays, and the quasinormal modes. 
To achieve this end we must have: (1) A perfect light 
absorber which can act as the event horizon and (2) A 
medium outside the black hole that replicates the black- 
hole spacetime. In the following sections we describe how 
a particular profile of permittivity and/or permeability 
of the optical analogue can lead to an optical spacetime 



that has the black-hole metric and therefore reproduces 
the gravitational phenomena outside the event horizon. 
Within the framework of our formalism, we further dis- 
cuss the broadband omnidirectional light absorber devel- 
oped by Narimanov and Khildishev [7| that serves exactly 
as a perfect absorber, or a black hole. We also show how 
one can, to a good approximation, simulate the black- 
hole spacetime with purely dielectric medium. 



II. ELECTROMAGNETIC WAVE EQUATION IN 
AN INHOMOGENEOUS MEDIUM: 

The macroscopic electromagnetic field strength tensor 
in a continuous medium is defined as[l[ 



(1) 



Where F 0jU = -F„ = = e ijk B k , P v = P 



LiU — ^ fi 7 X ij — MJA: J -- , j 1 — 1 

D v , P 3 = e ljk H k . E.ereF aT is the field strength tensor 
in free space and Z^ var is a fourth rank tensor whose 
components are the permittivity and permeability ten- 
sors. It is antisymmetric in the first two (/it, v) and the 
last two (tr, r) indices. e 1 ^ is the Levi-Civita symbol. The 
roman indices span 1, 2, 3 while the greek indices range 
over 0,1,2,3. The tensor Z is related to the permittivity 
and permeability tensors in the following way 



ryijkl -ijm kin.. — 1 



(3) 
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e y is the permittivity and fi mn is the permeability. Equa- 
tion (1) is the succinct form of the two known equations: 
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D = sE and H = yT^B . For a source-free medium in 
fiat space-time, the Maxwells equations involving sources 
are given by 



{Z»" rT F„\ v = iM i J' 1 
In the absence of sources 

(Z^ T F aT ) tU = 



(4) 



(5) 



For a curved space-time with metric g^ v , Maxwells equa- 
tions are given by - 



d 



(6) 



The similarity of eqs (5) and (6), suggest an analogy be- 
tween an inhomogeneous medium and a curved space- 
time, which is given by : 
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g^ a is a function of the permittivity and permeability and 
is used with a bar on the top to distinguish it from the 
actual metric. Here 7 is the determinant of the spatial 
metric in the ambient curvilinear co-ordinate systcm(r, 
0,4>) defined by j a p = -g a p+ (go a goi3)/goo , where g is the 
metric of the ambient flat spacetime. In spherical polar 
co-ordinate system on fiat space-time, ^/j — r 2 sind . 

One can directly formulate an optical spacetime at 
this stage itself, by defining an effective optical metric 
through the tensors e and fi , which are the components 
of Z. However, we proceed to derive the wave equation for 
a spherically symmetric, magnetisable dielectric medium. 
With this equation as the basis one can easily get the 
Hamilton- Jacobi equation for the geometrical optics limit 
and also analyse the case of purely dielectric medium. 
The formulation of an optical spacetime will be achieved 
through the similarity between the wave equation in an 
optical medium and that in a curved spacetime respec- 
tively, which will formally bring out the analogy between 
the metric and the permittivity/permeability. 

The field strength F ar is written in terms of the four- 
potentials , which are expanded in the vector spherical 
harmonics as follows M 



p — A _ A 
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A, = £K°' 0, ~a lm sinfl^) 



(8) 



smf 



00 



Odd parity funtions 
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Even parity functions. 

a lm , f lm , h lm , k lm are functions of r and t. Y lm { 



scalar spherical harmonics. Using the above in eq (5) 
gives, for odd parity functions the following equation 



i^ee a ,r ),r £ 



d 2 a lm 00 + 1) u 

Q t 2 r 2 



(10) 



And for even parity functions the following 



J)0 



d 2 b lm i{i + i) u i 



dt 2 



b lm = (11) 



Where b lm is given by e ee ^[{h lm ) fi - (f lm ), r ] = 

These equations are similar to those developed for the 
gravitational case, by J.A.Wheeler and D.Brill|9] and also 
by R.Ruffini, J.Tiomno, C.V.Vishveshwara[8J (who in- 
clude the source term) which are shown below 



{g rr af), r - g° 



d 2 a lm (1(1 + I)) j. 
dt 2 r'< 



a lm = 



(7) And for even parity functions it is 



(9 ), r -g 



no 



2ldm 



d 2 b 



1(1 + 1) 



dt 2 



b lm = 



(12) 



(13) 



These equations resemble the Regge- Wheeler equa- 
tions for gravitational perturbations. The similarity in 
equations (10), (11) and (12), (13) gives rise to similarity 
in a number of phenomena which depend on the form of 
the wave equation. One can formulate an optical space- 
time with an effective optical metric whose components 
are the permittivity and the permeability tensors. 

9 rr = Hsi = Hi 9 00 = -e 66 = S», = 1 



Metamaterials can be used to realize such an optical 
spacetime. Such a system has been studied particularly 
for the Schwarzschild spacetime in [l(j ■ 

Choosing e ee = (1 - ^y 1 and p,f£ = (1 - £),one gets 
an exact analogue of Schwarzschild black-hole spacetime, 
where r=L can be considered to be the location of the 
event horizon co-ordinate. The analysis of the gravita- 
tional bending of light and quasinormal modes will follow 
exactly as in the case of the gravitational case. 



III. WAVE EQUATION IN A 
NON-MAGNETISABLE DIELECTRIC MEDIUM: 

Now consider the case of a non-magnetisable dielectric 
medium (with /I = l).This is of interest because of the 
simplicity of using only a dielectric medium for experi- 
ments. 

With p, = 1 equation (10) becomes 



_ / Ira \ 
~=90\ a ,r ),r 
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da 



1(1+1) 



dt 2 



a = (J 



(14) 



To transform the equation to Regge- Wheeler form, sub- 



stitute a 




into the above equations to get 



P(r)}i/j = (15) 
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Here 



P(r) 
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So, the effective spacetime metric can be written as fol- 
lows 3 00 = , g rr = 

The function P(r) represents the effective potential, 
which in this case is not similar to that in (10). There- 
fore the behaviour of EM waves in the medium is slightly 
different but in certain approximation it will be similar 
to the gravitational case. The derivatives of e can be 
very small and hence we can neglect them in P(r). Further 
approximations depend on the specific form of e 66 con- 
sidered and also on the features to be studied using the 
analogue spacetime. 



IV. DISCUSSIONS: 

We discuss the application and extension of the formal- 
ism in the case of purely dielectric medium which will in- 
volve some approximations. For magnetisable dielectric 
medium the procedure is exactly the same. 



A. Schwarzschild Black-Hole analogue in dielectric 
medium 



Choice of e" 



gives us the analogue of the 



Schwarzschild spacetime. Substituting it into the wave 
equation gives 

9 \h Mi/ n qq + i)) ,., r *(* + i) r 



derivativesofe 60 ] = 



(16) 



Here L represents the Event Horizon of the analogue 
black hole. We are only interested in the physics of the 
phenomena happening outside the event horizon i.e. at 
r > L. If L is sufficiently small, then outside the black 
hole we can ignore terms like L/r 3 in the effective po- 
tential P(r).Also, if we consider a medium with its per- 
mittivity very slowly varying along r, we can ignore the 
terms with the derivatives of e 98 .Therefore we have a 
simplified equation 



d 2 jj 



(1 



L 1(1 + 1] 



(17) 



This is exactly similar to the gravitational case. Here 
we have used the Regge- Wheeler tortoise co-ordinate 
r* = r + L In ( — 1). The following are the two met- 
ric components of the effective optical spacetime g 00 = 
-(1 - ^y 1 , g rr = (1 - £). The above identification 
ensures that the optical spacetime in fact corresponds to 
the Schwarzschild with the event horizon located at r = 
L. Thus the optical system devoid of any magnetisable 
media can be used to simulate the black-hole spacetime, 
without the use of metamaterials. 



B. Geometrical optics approximation: Geodesies of 
light rays 

Once an effective optical spacetime is established, 
other phenomena like bending of light in the black-hole 
spacetime, must directly follow. Let us start with the 
cikonal substitution ijj — e s ^ r,t ' to the wavefunction in cq 
(15) .After some simplification we get 



P(r) = 



(18) 



Except for the effective potential term, this is similar to 
the Hamilton- Jacobi equation for light in curved space- 
time 







(19) 



Taking the specific case of the Schwarzschild metric, 
with the same approximations as discussed before we get 



1 r> y dt> K r' y dr 



(20) 



Here J is the angular momentum. The similarity in the 
Hamilton- Jacobi equations immediately implies similar 
geodesies. The equation of the trajectory in such a spher- 
ically symmetric spacetime is given by 



Jdr 



Co 



(21) 



In the effective optical medium the metric components 
are replaced by respective functions of permittivity, 
which act as effective metric. 



The optical Black Hole of Narimanov and 
Khildishev. 



In order to simulate black-hole phenomena, we must 
have (a) a medium which simulates the black-hole space- 
time and (b) a perfect absorber, which acts like a hori- 
zon, absorbing the entire incident light. We have al- 
ready developed the formalism for designing the equiva- 
lent optical spacetime around a black hole. The omnidi- 
rectional, broadband optical absorber of Narimanov and 
KildishevQ can be used to mimic the black hole. Their 
model consists of: 1. A central core of complex permittiv- 
ity, which absorbs all the incident light. 2. An outer shell 
, with some permittivity profile, which bends the incident 
light, guiding it to the core and thus ensuring that the 
entire incident light converges to the core, without any 
loss. The model with the core and the shell is placed in an 
external medium. Our formalism can be used to arrive 
at the Hamiltonian used by Narimanov and Kildishev 
to obtain a Keplerian power law variation of permittiv- 
ity which ensures proper convergence. Starting from our 
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cikonal equation (18) and assuming S = — Bi+function 
of r,0,cj> and ignoring the derivatives of permittivity we 
get 



i as 2 j 2 _ 2 

In this is expressed as 



1 J 2 



(22) 



(23) 



(Note that the authors in Q use "m" instead of "J" for 
the angular momentum.) This gives a ray trajectory 



, * (24) 



The permittivity choice e ee = e ee [l + (f )"] for n > 2, 
in the outer shell will enable maximum convergence and 
correspondingly facilitates maximum absorption. 

This concludes our discussion of the analog black-hole 
formalism and its application to spherically symmetric 
spacetimes. Our further investigations will include ex- 
tending the formalism to the developing of analogues of 
axially symmetric spacetimes and studying physical ef- 
fects in their background. 
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